Topological aberration of optical vortex beams and singularimetry of dielectric 
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The splitting of a high-order optical vortex into a constellation of unit vortices, upon total reflec- 
tion, is described and analyzed. The vortex constellation generalizes, in a local sense, the familiar 
longitudinal Goos-Hanchen and transverse Imbert-Federov shifts of the centroid of a reflected opti- 
cal beam. The centroid shift is related to the centre of the constellation, whose geometry otherwise 
depends on higher-order terms in an expansion of the reflection matrix. We present an approxi- 
mation of the field around the constellation of increasing order as an Appell sequence of complex 
polynomials whose roots are the vortices, and explain the results by an analogy with the theory of 
optical aberration. 

PACS numbers: 42.15.Fr, 42.25.Gy, 42.50.Tx 



Understanding and manipulating the spatial structure 
of light beams is a fundamental theme of modern optics. 
Beyond the simple structure of plane waves, beams with 
inhomogcneous complex amplitude can carry quantized 
orbital angular momentum (OAM) which has many 
applications, including free-space classical and quantum 
communication 043, superresolved microscopy [B| and 
extra-solar planet detection d, 0] 

Free space optical modes carrying quantized OAM 
have optical vortices (phase singularities) on their axis 
[1, Q: they have the amplitude structure r'*' exp(i£0) 
close to the beam axis, for t an integer: around the 
axis, as the intensity becomes zero, the phase gradient 
becomes infinitely large. The core region close to the 
axis of such beams is therefore exceptionally sensitive to 
any imperfections in optical components along the path 
of the beam. Any disruption to the pure azimuthal phase 
structure around the axis of the beam breaks the axial, 
stren gth- ^ vortex into a constellation of \l\ unit strength 
zeros [10( . We call this effect topological aberration, as the 
effect of aberration disrupts the topology of the simple 
optical mode. 

Here, we describe in detail the topological aberration 
which an OAM-carrying beam experiences under dielec- 
tric reflection by an oblique angle. We will see that the 
resulting constellation of zeros from the original strength- 
£ vortex depends on an aberration-like analysis of the 
complex reflection coefficient in Fourier space, which is 
sensitive to terms up to order \£\ in its Taylor expansion: 
it is therefore possible, using the spatial distribution of 
phase singularities, to determine higher orders of the ef- 
fective reflection-induced aberration of the beam. Fur- 
thermore, this phenomenon will affect any vortex beam 
reflected by a mirror, prism or beam-splitter; as measure- 
ment techniques become more sophisticated, apprecia- 
tion of this universal azimuthal symmetry breaking will 
become increasingly important. Although the fact that 
a high-order vortex breaks apart upon perturbation has 



been long understood, our main result here is a formula 
for aberrated constellation of vortices. 

Tracking the change of position of optical vortices 
on reflection echoes the study of optical beam shifts, 
where the centroid of a homogeneously-polarized beam 
is shifted by an amount proportional to the optical wave- 
length, either in the plane of incidence (Goos-Hanchen 
shift U|), transverse to it (Imbert-Federov shift [l2l.[l3|). 
or some combination of the two. The net shift requires 
the initial transverse polarization to be homogeneous, 
and that the beam be narrow and axisymmetric in inten- 
sity, but is otherwise independent of the original ampli- 
tude profile. In recent years, there has been an explosion 
of interest in various types of beam shift, including both 
fundamental insight and possible applications 14 1. 

It is simple to see the connection between beam 
shifts and aberration: for total reflection in the plane 
of incidence (p-polarization) or perpendicular to it (s- 



polarization) , the reflection coefficient [151 ] is a simple 
'tilt' [lj§ to first order; as the reflection coefficient acts 
on direction, in real space the reflected beam is trans- 
lated. In partial reflection when the modulus of the re- 
flection coefficients vary, there are also angular shifts of 
the beam's propagation 17 1. 

The notion of using the displacement of an on-axis zero 
to study beam shifts goes back to Hans Wolter in 1949 
[Tij wn0 observed that the position of a zero is a precise 
marker depending on local information, rather than the 
beam centroid, which depends on the global distribution 
of intensity. Our approach may be thought of as gen- 
eralizing Wolter's approach to any incident polarization 
and vortex order, all of which may be further extended to 
other aberrations of vortex beams, including refraction. 

Optical vortices are objects of scalar optics, and hence 
it is natural in our framework to consider a polarized 
component of the reflected beam; beam shifts for polar- 
ized components are analogous to quantum 'weak val- 
ues' of the reflection operator [lij], which have shifts in 
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FIG. 1. (Colour online) Incident and reflected beam geometry, 
a) Schematic of the relation between interface plane and our 
coordinate frame centred at the virtual reflected beam (red), 
with intensity profile given in b). The green cone depicts a 
section of constant d of the incident beam, the blue cone is 
the shifted, real reflected beam, with intensity as in c) for 
n = 2/3,8o = 44°. d) Schematic of the beam coordinates. 
The local plane and angle of incidence depend both on 5 and 
a. e) Beam coordinates for a I = —4 vortex beam (phase 
indicated by colour wheel) and local resolution into e s and 



both real and Fourier space. The spatial shift to the 
ccntroid of vortex-carrying beams involves both 'spatial' 
and 'angular' parts of phase-flat beams wc will 

see that this general shift is simply the first term, analo- 
gous to tilt, in an aberration expansion of the polarized 
component's reflection matrix. In our calculation, higher 
aberration orders of the reflection matrix are realized as 
coefficients of a complex polynomial approximating the 
low amplitude in the core of the reflected beam: the vor- 
tex positions are related to these aberration terms by the 
fundamental theorem of algebra. 

The basic beam geometry is demonstrated in Fig. [T] 
choosing a natural coordinate system based on the con- 
struct of a virtual reflected beam 23[ , which is obtained 
by specular reflection of every plane wave in the incident 
beam. The central propagation direction is (0,0, 1), and 



we will make much of spherical angles about this direc- 
tion: azimuth a and polar angle S. Since reflection flips 
the sign of £, the reflected beam is a(S) exp(— i£a), with 
real a(5) tightly centred around 5 = 0, such that the sec- 
ond moment (72 = J dS\a(5)\ 2 S 3 <C 1 (where the extra 
5 ss sin 6 is the jacobian in spherical polars). The upper 
limit of this integral is larger than the width of the spec- 
trum. In these coordinates, the normal to the incidence 
plane has direction (— sin 9q, 0, cos 9q), so, for each com- 
ponent of the spectrum labelled by 5, a, the incidence 
angle 9 is given by cos 9 = cos 5 cos 9q ~ cos a sin 8 sin 9q . 
This angular dependence appears in the reflection coeffi- 
cients r s (6) and r p (9) [l6| , which are defined with respect 
to the local plane of incidence for every plane wave within 
the spectrum (see Fig. QJi). The indices s and p thereby 
refer to the directions orthogonal (s) and parallel (p) to 
the local plane of incidence (see Fig. QJ) and the general 



reflection matrix R 



P p consists of the projec- 



tors Pj = ej ej,j = s,p of the incident field onto the 
local e s and e p direction (see Fig. QJ) multiplied by the 
appropriate reflection coefficient. R acts on an initial po- 
larization E, which, for the central wave vector, is a 2D 
constant Jones vector with components E x and E y . This 
is the global polarization of the incident field which we 
distinguish from the local polarization in terms of s and 
p components. Because of transversality, however, the 
other plane waves in the spectrum have a small z compo- 
nent such that E = (E x ,E y ,— [E x cos a + E y sin a] tan <5) . 
As we only consider a polarized component of the beam, 
the reflected beam is filtered by a constant polariza- 
tion analyser F with components (F x ,F y ,0). All of the 
physics of beam reflections can be explained by a Taylor 
expansion about 5 = of the scalar multiplication oper- 
ator F* ■ R ■ E in Fourier space as we briefly summarize 
in the following (also see 24[). 

In the transverse plane the real-space scalar reflected 
beam is given by 



ip(r) = dS daa(6)6e lk 



ik sin 5r- (cos a, sin a) — i£c 



F* R E, 

. (1) 

where r = {x, y) and k is the wavenumber of the incident 
light. The formula for the beam shift follows from a first 
order Taylor expansion of the reflection matrix term, 



F* R E ps R[l + 5T> ■ (cos a, sin a)] 



J^q3T>- (cos a, sin a) 

(2) 

where D = {T> x ,T> y ) is a complex 2- vector with compo- 

-F;E y r s ) and 
-r s )/{F*E x r p -F;E y r s ). 



nents V x = (F*E x r' p - F* y E y r' s ) / \F*E x r p - 
V y = {F;E x + F*E y )cote (r 

When £ = 0, the spatial shift is given by -ImD/fc, 
which adds to the {x, y)-dependent term in the exponent 
in (JTJ when identifying sin S with S. This is the well- 



known Artmann formula for spatial beam shifts 25| . The 
angular shift is the angular mean of the Fourier transform 
of ip, which is given by o^Rel? 0,[2^|. For £ ^ the net 
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FIG. 2. (Colour online) a) Plot comparing the constellation 
of vortices obtained from the roots of the local expansion in 
pjl with a numerical calculation of the phase of the field. 
The incident field is a Bessel beam incident at Go = 46° with 
I — —4 with a fixed opening angle of 5 — 0.01. The incident 
polarization E and the analyser F are both oriented in the 
x direction (TM/TM). b) Plot showing the variation of con- 
stellations for 45° diagonal incident polarisation and different 
analyser settings including linear in the y direction (TE) , right 
(RC) and left (LC) circular polarization, as well as 45° and 
135° diagonal polarization. 



vortex beam shift is 



D 



ccntroid 



-ImX> - \£er 2 ■ Rc2>, 



(3) 



where er 2 is the second Pauli matrix 



llliil. The first 
term here is the usual Artmann translation; the second 
term, usually associated with the angular shift, comes 
from the azimuthal complex amplitude structure of the 
vortex beam. 

The shift to the intensity centroid of a vortex-carr yin g 
beam is different from the shift of the vortex itself [21. 27|. 
The translation — ImZ> affects the entire beam including 
vortex and the centroid. However, ReX>, being associated 
with a change in the profile of the beam, affects them 
differently: for an asymmetric profile, a vortex, being an 
absence of intensity, repels the intensity centroid, and it 
follows form the more general argument below that the 
shift of a vortex in a \t\ = 1 beam is 



ImX> ± ier 2 • ReX>, 



(4) 



where ± refers to £ ^ 0. 

When \l\ > 1, topology preserves the overall vortex 
charge, but the symmetry is broken and typically there 
is a constellation of \£\ unit-charge vortices in the re- 
flected beam, with separation comparable to the beam 
shift itself. An example of a numerical calculation of the 
reflection of a £ = —4 Bessel beam (28|, with several dif- 
ferent incident and analyzer polarizations, is shown in 
Fig. [21 the constellation is not simply a regular polygon 
or 'row' [111, but a more complicated function of incident 
and analyzer polarization, incidence angle 9q and refrac- 
tive index n. The ccntroid of the vortices, represented by 



the white circle, is the same as the position of the sin- 
gle shifted vortex D vor te X : the centroid of the \£\ vortex 
points is a discrete, topological counterpart to the shift 
of the intensit y c entroid, but without the £ weighting, 



anticipated in [21 



To analyse these constellations we now derive an an- 
alytic approximation for the reflected vortex beam close 
to the beam axis, as a complex polynomial in x + iy or 
x — iy depending on the sign of £ [29|, by collecting all 
the terms of the lowest order in 6 in the expansion of the 
reflection matrix ©. Up to order \£\, this can be written 
as 



F* R E « R I 1 + SCi + ls 2 C 2 + 



1 



(5) 



where R = F* ■ R ■ E\s=o- The original beam shift fol- 
lows from C\ = T> ■ (cos a, sin a), which contains only 
single powers of sin a and cos a and depends on the first 
derivatives of the reflection coefficients. The higher co- 
efficients C n contain combinations of sin u a cos v a and 
— (d q / d8 q )r j\$ = Q for j = s,p with u + v = n and 
q < u,v. Each C n may thus be written as a complex 
Fourier series 



c„ e 



2" 



(0) 



This expansion of the reflection matrix is analogous to an 
aberration expansion in terms of complex Zcrnikc modes 
[30| . with n giving the order of aberration: corre- 
sponds to tilt, C2 to astigmatism, etc. 

The core of the reflected vortex beam can be found by 
combining Eqs. ([T]) and (0; after integration over the az- 
imuthal a, each c ±lqa contributes a Bessel function Je^q, 
which is Taylor approximated close to the axis. The low- 
est order term in the beam is <5' £ ', which comes from 
each for £ ^ 0; all other terms have higher order and 
therefore smaller contribution. The size of the nth con- 
tribution, corresponding to is 



z±5 n 



2"n! 



dac ik5r cos ( Q - ( A)- i (^ = F™) Q 



5 w c±2niW- n [2^n\{\£\ - n)!]" 1 [k(x =F iy)] 



(7) 



substituting r exp(=Fi</>) = x, =fiy. 

On collecting common factors, the reflected field near 
the axis for a given £ is proportional to 

^ocC |f| -iK|c^ | - 1 -^^c 2 t C l ^ 2 + --- + (-i) W ^ | , 

(8) 

where we have used ( = k(x =F w) with '— ' for £ > 
and '+' for £ < 0. Each root of this complex polynomial 
corresponds to a vortex, and the constellation depends 
on the roots of the polynomial with coefficients given 
by (^i)" ™ times a binomial coefficient as in ((SJ). For 



4 



increasing 



l\ the polynomials in ([5J form an Appell se- 
quence [3J|, with coefficients generated by the sequence 
c n of coefficients from the Zernike expansion of the reflec- 
tion matrix, for which there are many relations between 
the roots and the coefficients [32[ . The integral over 6 as 
in dU) docs not affect Eq. ©, since these terms have the 
same, lowest-order contribution in 5, and so the constella- 
tion is independent of the spectrum and hence the radial 
profile of the beam, as long as it is narrowly confined in 
Fourier space. 

In the complex plane, the arithmetic mean (centroid) 
of the roots of a polynomial of the form (|8j) is given by icf 
[3^ | , which proves that the mean position of the vortices 
is always given by icf = "D ■ (i, ±1). The vortex centroid 
shift is independent of £, and can be understood purely in 
terms of the representation of the beam as a polynomial 
in the complex (x =p iy) plane. 

The vortex constellation for \£\ > 1 is sensitive to terms 
beyond the first-order shift. The second order term cf 
is related to the complex astigmatism of F* HE, pro- 
portional to e ±2iQ . When \£\ = 2, Eq. (|8]) is a quadratic 
equation with the 'constant'; icf is the midpoint of 
the two roots, and their separation from the midpoint is 

2 — (cf ) 2 (where the complex argument represents 
the direction) . Measurement of the two singularity posi- 
tions would therefore provide the two complex numbers 
cf and C2 , from the first two orders of the Zernike ex- 
pansion of the reflection coefficient. The former contains 
zcroth and first derivatives of the reflection coefficients 
via the complex shift vector T> and the latter is given by 



± 



c? = (F*,F;)-cf-(E x ,E y )/(F* 
the C2 matrix is given by 



E x Tp 



F*E y r s ), where 
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FIG. 3. (Colour online) Vortex constellations for increasing 
positive I (a) and negative I (b). c) Contours in Fourier space 
(conical beam geometry) of the real and imaginary part of 
F* • R ■ E for 6» = 46° and both F and E set to 45° diagonal 
polarization, d-g) Contours plots for of the real parts of the 
coefficients in © multiplied by the appropriate Fourier factor 
exp(— ina) for comparison with I < in b). d) n = 1, e) 
n = 2, f) n = 3, g) n = 4. (Scaling identical to figure c). The 
gradient over the contours in d) corresponds to the shift as 
indicated by the black line between x and o in b). 



(r s +r p ) cot 2 6*o 



I ±i 

±isc 2 -1 



Ti(r' s +r' p )cote 



I 

1 



■ cot ( 



Zl'P 



(9) 



where sc 2 = sec 2 9<j. This complicated form of shows 
that for higher order shift effects a separation in terms of 
diffractive corrections and optical spin-orbit interaction 
is no longer possible 
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In Fig. [3] we illustrate how this 'singularimetry' would 
give increasing information as positive I increases (a) and 
as negative i increases (b). As \£ | increases, the constel- 
lations increase in complexity, but maintaining common 
features (such as their centre), as their coefficients are al- 
gebraically related in Eq. . These complex coefficients 
are proportional to the c„ , which correspond to succesive 
azimuthal aberration terms (tilt, astigmatism, ...). This 
is shown in Fig. [3] c) - g), which illustrates how a fixed 
reflection matrix term F* ■ R • E as a function of a, 6 (c) 
is decomposed into aberration terms of increasing order. 

It has recently been appreciated that there is a close 
connection between beam shifts and the notion of quan- 



tum weak values [33(. We note that going beyond \£\ = 1 
shows that there is weak information beyond the first 
order weak value. 

Our argument extends beyond reflection, to any parax- 
ial vortex beam experiencing an aberration-like multi- 
plication in Fourier space: rather than a single order I 
on-axis vortex, the beam will have a constellation of I 
vortices which can be represented by a polynomial with 
coefficients dependent on the azimuthal terms. This ef- 
fect of topological aberration will affect any vortex beam 
undergoing oblique reflection, refraction in a real opti- 
cal system, suggesting a fundamental limitation on the 
purity of quantized OAM modes in optical devices 34 [ . 
However, as the connection between vortex positions and 
the aberrations is so algebraically direct, any singulari- 
metric device which measures the precise vortex distribu- 
tion has direct access to the beam shift and later terms, 
without any further analysis of the overall beam. 
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